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We investigate the three-dimensional 0(2) model on lattices of size 8'^ to 160'^ close to the critical point at 
zero magnetic field. We confirm explicitly the value of the critical coupling Jc found by Ballesteros et al. and 
estimate there the universal values of Qr and £,/ L. At the critical point we study the finite size dependencies of 
the energy density e and the specific heat C. We find that the nonsingular part of the specific heat Cns is linearly 
dependent on 1/a. From the critical behaviour of the specific heat for T ^ Tc on the largest lattices we determine 
the universal amplitude ratio A'^ /A~ . The a-dependence of this ratio is close to the phenomenological relation 
A+/A- = f - 4a. 
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1. INTRODUCTION 

0{N) models in three dimensions play an im- 
portant part in condensed matter physics and in 
quantum field theory, because many systems be- 
long to the corresponding universality classes. In 
three dimensions the case = 2 is a special 
one, because the 0(2) model is the first 0{N) 
model (with increasing N) exhibiting massless 
Goldstone modes. Furthermore its critical expo- 
nent a, which controls the critical behaviour of 
the specific heat, is negative and very close to 
zero. In the famous shuttle experiment the 
universal ratio A'^ /A^ and a have been deter- 
mined experimentally. Here we want to calculate 
this ratio from Monte Carlo simulations. 

The 0(2)-invariant nonlinear tr-model (or XY 
model) for zero magnetic field, which we examine 
here, is defined by the partition function 



[d$]exp[J 
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Here c/Jx is a 2-component unit vector at site x 
and J = l/T is the inverse temperature. We use 
the lattice average of the spins 
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to define the magnetization M, the order param- 
eter of the system, by 



M = (|rn|), 

the susceptibility x 

and the Binder cumulant 
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The second moment correlation length £^2nd is cal- 
culated from 



2nd 
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where F is the Fourier transform of the correla- 
tion function at momentum = 27re^/L. Im- 
portant observables for this work are the energy 
E, the energy density e 
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and the specific heat C 
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At the critical coupling Jc the finite size be- 
haviour of the energy density is 



e{L) = ens+QcL " , 

and that of the specific heat 



(9) 
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ens and Cns are the nonsingular parts of the en- 
ergy density and of the specific heat, respectively. 
In the thermodynamic limit the critical behaviour 
of C for T close to is 

cit) = Cns + — ir" [1 + 4 \tr + 4 1], (11) 

a 

where t = {T — Tc)/Tc is the reduced tempera- 
ture. Here correction to scaling terms have been 
included. We use the first irrelevant exponent 
Lo = 0.79(2) from ||^ in the following fits. 

2. NUMERICAL RESULTS 

Our simulations were done on three-dimensio- 
nal lattices with periodic boundary conditions 
and linear extensions L = 8 — 160 using the Wolff 
single-cluster algorithm. 

First we have determined again the critical cou- 
pling Jc, utilizing the fact that Binder's cumu- 
lant should be finite size independent at critical- 
ity. We have interpolated data from simulations 
on lattices with L = 24,36,48,72 and 96. The 
curves for the different lattices do not cross in 
a single point due to small corrections to scaling. 
Including these we find that the shift A from crit- 
icality of the crossing point between two lattices 
of size L and L' {h ~ L' / L) is 



AJ, 



L.L' 



j^ — LO — l/ U 
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The i^-dependence of AJ^'-^ is not relevant as 
long as S [0.669; 0.675]. Extrapolating to L ^ 
oo we find the value 



= 0.454167(4), 



(13) 



in agreement with Jc — 0.454165(4) from In 
a similar way we determine the universal values 

of Qr and S,2nd/L at Jc to 

= -1.758(2) and 6nd/i = 0.593(2), (14) 



in accordance with |Q] and . 

In order to estimate the nonsingular parts ens 
and Cns of the energy density and the specific 
heat the finite size effects of these observables at 
Tc were studied. The model was simulated at 
Jc = 0.454165 on lattices with L = 8 to L = 160. 
Fits to the data with Eqs. (^ and (^ show no 
corrections to scaling in case of the energy den- 
sity and only small corrections for the specific 
heat. The quantity e„s exhibits no noticeable de- 
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Figure 1. The energy density (a) and the specific 
heat (b) vs. L at criticality. The dashed line 
shows e„s and the solid lines fits for i' = 0.671. 

pendency on h', and we find e„s ~ —0.98841(3). 
When we treat as a free fit parameter, we get 
0.671(2). In case of the specific heat the situa- 
tion is different. Its nonsingular part varies from 
about 50 for v = 0.669 to 16 at = 0.675. This 
is so because the exponent a/iy in Eq. (|lO| ) is ap- 
proximately zero. As shown in Fig. H the non- 
singular part of the specific heat Cns is linearly 
dependent on 1/a. We find 



Cnsia) = 3.35(26) - 



0.3176(43) 



(15) 
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Figure 2. The nonsingular part of the specific 
heat. The dashed Hne is the fit ([l5|). 

This result for C„s is used to analyse the specific 
heat at T Tc and ^ oo. Our data are shown 
in Fig. ^. We have interpolated the data from 
the largest lattices by reweighting. The resulting 
curves have been fitted to the form (|ll|). In the 
broken phase the correction to scaling terms in 
the bracket are relevant for the fit, whereas these 
terms are negligible for T > Tc- Again we observe 
an a-dependency in the fit parameters, especially 
for the amplitudes A^: 

A+ = 0.3177(2) - 3.29(4)a+ 18. 8(1.3)a2(ig) 
A- = 0.3176(3) - 1.97(4)a + 7.8(1.4)a^ (17) 

As expected, we obtain the same amplitudes for 
a ^ 0, and the 1/a-pole term in C„s is cancelled 
exactly there: the same specific heat data can as 
well be fitted at a = 0. 

The universal ratio A'^ /A^ can now be written 

as 



A+/A- = 1 -4.23(3)a-f 3.3(1. 8)a2 



(18) 



It is shown as solid line in Fig. ]^ This result is 
well in accordance with former results, e.g. from 
the shuttle experiment [|l| and analytic determi- 
nations Ql and The leading part is also close 
to the phenomcnological relation [|j 
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Figure 3. The specific heat for different L vs. the 
coupling J. The line shows the position of Jc- 
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Figure 4. The universal ratio jA as a func- 
tion of —a. The solid line shows the result (p^. 
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